Measurements of Grain Motion in a Dense, Three-Dimensional Granular Fluid 



Xiaoyu Yang, Chao Huan, and D. Candela 
Physics Department, University of Massachusetts, Amherst, Massachusetts 01003 

R. W. Mair and R. L. Walsworth 
Harvard-Smithsonian Center for Astrophysics, Cambridge, Massachusetts 02138 

(November 29, 2001) 



We have used an NMR technique to measure the short- 
time, three-dimensional displacement of grains in a system of 
mustard seeds vibrated vertically at 15g. The technique aver- 
ages over a time interval in which the grains move ballistically, 
giving a direct measurement of the granular temperature pro- 
file. The dense, lower portion of the sample is well described 
by a recent hydrodynamic theory for inelastic hard spheres. 
Near the free upper surface the mean free path is longer than 
the particle diameter and the hydrodynamic description fails. 

PACS numbers: 45.70.-n, 81.05.Rm,76.60.Pc 

Under suitable conditions, a noncohesive granular sys- 
tem such as dry sand can be maintained in a fluid-like 
state by continuous vibration of the container. A ques- 
tion of current interest is the extent to which such a sys- 
tem can be quantitatively described by a hydrodynamic 
theory, generalizing conventional nonequilibrium statisti- 
cal mechanics for elastic systems such as liquids and gases 
■ The key distinction of granular systems is the con- 
tinuous loss of kinetic energy during inelastic collisions of 
the grains, resulting in highly out-of-equilibrium states 
with strong gradients or time derivatives of macroscopic 
variables such as density and granular temperature. 

Theoretical and experimental studies of the "vi- 

brofluidized" state usually have treated simplified con- 
ditions of uniform spherical particles with nearly elastic 
collisions, often in reduced dimensions. It is important 
to determine the extent to which theories developed for 
this idealized domain can describe more typical, techno- 
logically relevant systems. Such systems often have ir- 
regularly shaped grains with a distribution of sizes, and 
complex energy loss mechanisms only loosely modeled by 
the usual restitution coefficient |l^ . 

In this paper, we present an experimental study of a 
three-dimensional system of mustard seeds fluidized by 
vertical vibrations of the container. Pulsed field gradi- 
ent nuclear magnetic resonance (PFG-NMR) combined 
with magnetic resonance imaging (MRI) was used to 
measure the distribution of positions and short-time dis- 
placements of the grains. In this way, we were able 
to measure both the density and granular temperature 
profiles of the system. To our knowledge, the only 
previous measurements of granular temperature profiles 
in three-dimensional vibrofluidized systems were the re- 



cent positron-emission particle-tracking (PEPT) studies 
of Ref. 0. The PEPT and NMR techniques are com- 
plementary as the former traces the motion of a single 
tagged particle over long times, while the latter measures 
nearly instantaneous averages of the motion of many par- 
ticles. The NMR method used here achieves better time 
resolution than previous PEPT work |l|] and is used 
to study a much denser system. Diffusing-wave spec- 
troscopy has also been used to study grain motion in 
three-dimensional systems l]!^ . 

Two-dimensional systems have been extensively stud- 
ied using video methods ||7|"pd||, and spin-tagged MRI 
has been used to measure much slower granular flows 
than studied here ||l^. Seymour, et al. have used PFG- 
NMR/MRI to image the diffusion of grains in a rotating 
drum and infer the granular temperature Jl6t . In contrast 
to that work we study a simpler physical system for which 
first-principles hydrodynamic theories are currently avail- 
able, and use more rapid gradient pulses to probe deep 
within the ballistic regime. In earlier work [ p^t we re- 
ported MRI measurements of the density of a vibrofiu- 
idized granular system. The important advance reported 
here is the ability to measure the granular temperature 
profile. Despite the small size and high density of our 
system (roughly 50% of the close-packing density), we 
find the density and granular temperature profiles to be 
surprisingly well described by a recently developed sta- 
tistical theory |Q for inelastic hard spheres (IHS). Near 
the top of the system, which is a free boundary, the mean 
free path becomes long and the hydrodynamic descrip- 
tion based on the IHS model breaks down. The horizontal 
and vertical degrees of freedom become decoupled, and 
the horizontal granular temperature grows with height 
despite "heating" only from below. We do not yet have 
an explanation for this phenomenon. 

The experimental setup was as follows. The sample 
cell was a vertical glass tube of inside diameter 9.00 mm, 
with a perforated Teflon bottom wall. It was sufficiently 
tall to prevent collisions of the grains with the top. The 
cell was supported at the center of the NMR probe and 
magnet (static field 0.975 T) and vibrated vertically by 
the arrangement described in Ref |l^ . The NMR proto- 
col used a stimulated-echo PEG sequence with a single 
gradient pulse pair followed by frequency-encoded one- 
dimensional MRI [Q. This allowed us to measure the 
joint distribution p(Au,w; At,t) of grain displacements 



1 




FIG. 1. Density n{z,t) of grains in arbitrary units as a 
function of height z and time t in the vibration cycle. The 
zero of the z axis is arbitrary. The thick sinusoidal line shows 
the amplitude and phase of the cell motion measured by the 
accelerometer. It is drawn at the position of the cell bottom 
as inferred from the density data. For most of the vibration 
cycle, the sample "floats" well above the cell bottom. 
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FIG. 2. Root mean square horizontal and vertical dis- 
placements Axrma (solid squarcs) and Azrms (open circles) 
averaged over the sample volume and time, as functions of 
the displacement observation interval At. The solid lines have 
slopes of 1.0 and 0.5, corresponding to ballistic and diffusive 
dependence respectively of the RMS displacement on At. 



Au along one axis during At and grain positions v along 
a second axis, as a function of At and the time t within 
the vibration cycle. By using very brief (150 /is) gradi- 
ent pulses, the displacement interval At could be made 
as short as 850 /is. The PFG sequence was not velocity 
compensated, so both collective motion (such as convec- 
tion) and random motion could be observed ]l6| , p^ . 

A complete four-dimensional data set p{Au,v; At,t) 
is large, and in addition the physical conditions such as 
number and size of grains, vibration amplitude and fre- 
quency, and ambient gas pressure can be varied. For this 
study we focussed on a single sample state, consisting of 
60 mustard seeds vibrated sinusoidally at 50 Hz and ac- 
celeration amplitude (14.85±0.15)g where g = 9.80 m/s^. 
The ambient pressure was reduced below 200 niTorr [ [l9| . 
The mustard seeds were roughly ellipsoidal with a mean 
diameter a — 1.84 mm, standard deviation of the diam- 
eter Act = 0.09 mm and mean mass m = 4.98 mg. 

We measured the joint distribution of horizontal or ver- 
tical displacements and vertical positions p{Ax, z; At, t), 
p{Az, z; At, t) as well as the transverse mass profile p{x) 
{x is horizontal measured from the cell center and z 
is vertical). To within the noise p{x) fits the form 
[1 — (x/R)^]^/^ for uniform transverse density in a cylin- 
der of radius R, showing that there was no large variation 
of particle density with radial coordinate. 

Figure ^ shows the vertical density profile as a function 
of time in the vibration cycle. The profile is largely in- 
dependent of t except for a disturbance to the lower part 
of the sample near the peak of the vibration waveform. 



when the cell impacts the sample. These time-resolved 
data suggest it might be possible to describe the bulk of 
the sample as a steady-state statistical system "heated" 
by energy input from the vibrating cell bottom. 

The RMS grain displacements Axrms , Azrms were 
computed from the decay of the signal magnitude vs. 
QxtQz, which are variables Fourier-conjugate to Ax, Az 
|p^ , p^ . Figure ^ shows these mean-square displacements 
averaged over the sample volume and observation time 
t, as functions of the displacement interval At. These 
graphs show a crossover from ballistic (Axrms, Azrms 
At) to diffusive {Axrms, Azrms oc Ai^/^)) time depen- 
dence. The diffusive power law is not clearly established 
by the data, but it is the ballistic-regime information that 
is used for the analysis below. 

The distribution of horizontal displacements p{Ax, z) 
is shown in Fig. ^, for At — 1.35 ms well within the bal- 
listic regime. Thus p{Ax/At,z) should approximate the 
distribution of horizontal velocities p(vx, z) It can be 
seen that p{Ax, z) is symmetric in Ax and roughly gaus- 
sian at each height z, at least well away from the sample 
bottom. This suggests that the motion is dominantly 
statistical with minimal large-scale convection. The dis- 
tribution of vertical displacements (not shown here) is 
nearly symmetric at the sample top, but becomes skewed 
as the sample bottom is approached. This reflects the 
asymmetry between upward moving grains that have re- 
cently impacted the cell bottom, and downward moving 
grains that have not done so. 

We fit the distribution of horizontal displacements 
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FIG. 3. Distribution p{Ax, z) of horizontal displacements 
Aa; and heights z of grains averaged over the vibration pe- 
riod, for displacement measurement time At = 1.35 ms. The 
vertical axis is arbitrary. 

p{/S.x,z) to the form exp — [(A2;/s2;)^/2]^ with Sx and (3 
as fitting parameters (/? = 1 gives a gaussian). The best 
fit (3 varies from 0.9 at the sample bottom to nearly 1.0 
at the sample top, although the deviation of (3 from unity 
is only marginally significant even at the sample bottom. 
Therefore the deviation of the velocity distribution from 
gaussian in this three-dimensional system is less than 
seen experimentally in two dimensions, for which one re- 
cent study found /3 « 0.75 The RMS grain displace- 
ments computed from the fits agree to within 10% with 
Asrms, Azrms computcd directly as described above. 

To compare the data with theory, the density was nor- 
malized using J" n{z)dz = ua where ua = 119.6 cm~^ 
is the areal density pO| . The granular temperatures for 
horizontal and vertical motion were computed from the 
ballistic-regime (At — 1.35 ms) mean-square displace- 
ments as Tx = m{Ax^)/At^, = m{Az^) / At'^ . Dimen- 
sionless height, number density and temperature vari- 
ables were computed from z* — z/a, n* = na^ and 
T* ^ — Tx.z/rnga using the values of a, m and g given 
above. Figure^ shows the measured n*(z*) and T* .,(z*) 
as well as the dimensionless mean free path l*{z*) = ija, 
estimated in a manner similar to Ref. pi]] . 

The maximum observed density n* k, 0.7 should be 
compared with the values n* = 1.2, 1.41 for poured ran- 
dom and crystalline close packed monodisperse spheres, 
respectively 22 1. In the lower part of the sample [z* < 



4) , the horizontal and vertical granular temperatures are 
nearly equal to one another and decrease monotonically 
with height. Although energy is supplied to the vertical 
degrees of freedom by collisions with the cell bottom, the 
mean free path is much less than the particle diameter 
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FIG. 4. (Top) Measured dimensionless number density n* 
(filled squares). (Center) Measured dimensionless horizontal 
granular temperature (filled squares) and vertical granular 
temperature (open circles) . (Bottom) Dimensionless mean 
free path £* , computed from the density data. These quan- 
tities are plotted versus vertical position measured in grain 
diameters, z* — z/a. The smooth curves show the number 
density n* and granular temperature T* calculated from the 
theory of Ref. [4], with the heat current at sample bottom 
9z(0) ~ 3.4 and restitution coefficient a — 0.90 adjusted to 
fit the data. The vertical dotted lines in the top graph show 
the extent of the granular sample in the absence of vibration. 



{£* ^ 1) in this dense portion of the sample, result- 
ing in effective equilibration between T* and T*. Above 
z* = 5 the vertical temperature T* continues to decrease 
slowly with height while the horizontal temperature T* 
increases rapidly, falling far out of equilibrium with T* . 
It is reasonable that T* and T* should become unequal 
where the mean free path is long, but we have not yet 
identified the physical mechanism that causes T* to in- 
crease rapidly with height as measured. 

The smooth curves in Fig. ^show a fit to the data of the 
hydrodyamic theory for inelastic hard spheres of Garzo 
and Dufty [Q , which is applicable for wider ranges of den- 
sity and inelasticity than other available theories. In in- 
elastic fluids the heat current is q = — kVT — /iVn where 
K is the thermal conductivity and is a new transport 
coefficient that vanishes for elastic systems [ p3||24| . We 
have numerically integrated the expressions for k, fj,, and 
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the pressure p from Ref. M with respect to z* to obtain 
the curves shown in Fig. m. Due to the open top bound- 
ary the dimensionless pressure at the bottom is p*(0) = 
{a'^ / mg)p{Q) = n\. Here n\ = a'^UA = 4.04 is the 
dimensionless areal number density. The upward heat 
current at the cell bottom g*(0) (CT^/^/g^/^m)gz(0) 
and the restitution coefficient a are used as fitting pa- 
rameters. Presumably q*(0) is set by the amplitude and 
frequency of cell vibrations. The density at the bottom 
n*(0) is set to give vanishing heat current at the cell top 
(q*(z*) ^ as z* +00). 

Once 52(0), a and the sample bottom location are ad- 
justed to fit the experimental density profile n*(z*), the 
granular temperature profile T*{z*) is predicted by the 
theory with no further adjustable parameters. As shown 
in Fig. ^ the vertical temperature T* agrees well with 
the theory over the entire measured range of z* . How- 
ever, the horizontal temperature T* deviates from the fit 
for z* > 5. We find no choice of fitting parameters that 
duplicates the rapid increase of T* seen experimentally. 
Therefore our data do not prove fi ^ 0, which is what 
allows the theoretical profile to have dT* /dz* > 0. 

Are the fitted parameter values g^(0) = 3.4, a — 0.90 
physically reasonable? For an open vibrated system the 
power input is Q = 2Nmg(V) where Nmg is the total 
sample weight and (V) is the average upward cell velocity 
during grain impacts . Combined with the fitted (0) 
this gives (V) = 5.65 cm/s = 0.12wo, where vq is the peak 
cell speed. This value for (V) appears reasonable as most 
of the collisions with the cell bottom occur near the top 
of the vibration stroke (Fig. |l|). 

We have attempted to measure a for the grains used in 
this study by videotaping their trajectories when dropped 
onto a hard surface. The apparent a varies from trial 
to trial, as the irregular shape causes conversion of some 
translational energy into rotational energy at impact. For 
an impact velocity of 44 cm/s (larger than is typical in 
the fluidized sample) the mean apparent a measured in 
this way was 0.79 and the largest value measured across 
all trials was 0.92. While not an accurate measurement 
of the restitution coefficient in the fluidized state, these 
results are at least consistent with the fitted value. 

In summary, we have used NMR methods to measure 
three-dimensional grain motion and density and granu- 
lar temperature profiles for a vibrofluidized granular sys- 
tem composed of irregularly shaped mustard seeds. The 
dense lower region of the sample is well described by a hy- 
drodynamic theory for inelastic hard spheres. Thus, the 
continuum theory appears accurate in spite of the small 
number of grains in the sample and the steep granular 
temperature gradient. In the upper region of the sample 
the mean free path becomes long, and the horizontal and 
vertical grain velocities become decoupled in a manner 
that is not described by the theory. 
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